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Abstract 

It is well known that the perturbation equations of massless fields for the Kerr-de 
Sitter geometry can be written in the form of separable equations. The equations 
have five definite singularities so that the analysis has been expected to be difficult. 
We show that these equations can be transformed to Heun's equations, for which 
we are able to use known technique for the analysis of the solutions. We repro- 
duce results known previously for the Kerr geometry and de Sitter geometry in the 
confluent limits of the Heun's functions. Our analysis applies can be extended to 
Kerr-Newman-de Sitter geometry for massless fields with spin and ^. 
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1 Introduction 



One of the most non-trivial aspects of the perturbation equations for Kerr geometry is the 
separabihty of the radial and the angular parts. It was Carter [1] who first found that the 
scalar wave function is separable in the Kerr-Newman-de Sitter geometries. Later, this 
observation has been extended for spin 1/2, electromagnetic fields, gravitational perturba- 
tions and gravitino for the Kerr geometries and even for the Kerr-de Sitter class of geome- 
tries. The equations are called Teukolsky equations[2]. Except for the electromagnetic 
and gravitational perturbations, the separability persists even for the Kerr-Newman-de 
Sitter solutions. An important application of this fact was the proof of the stability of 
the Kerr black hole [3]. 

Though the Teukolsky equation for Kerr geometries is separable, both angular and ra- 
dial equations has two regular and one irregular singularities so that the solutions cannot 
be written by a single form of any special functions, but they are expressed as series of 
special functions whose coefficients satisfy the three term recurrence relations [4- 7]. The 
solution of the angular equation is expressed in the form of series of Jacobi functions [4]. 
The solution of the radial equation is rather complicated because we need the solution 
which is valid in the entire region extending from the outer horizon to infinity. The 
solution is written in the form of series of confluent hypergeometric functions [5] which 
is convergent around infinity and the hypergeometric functions [6- 7] which is convergent 
around the outer horizon. By matching these solutions in the region where both solu- 
tions are convergent, we can obtain the solution which is valid from the outer horizon to 
infinity[6-8]. The great benifit of this kind of solutions is that the coefficients of the series 
are obtained as the Post Minkowskian expansion [6- 7]. This technique has been success- 
fully applied for the Post Newtonian expansion of gravitational waves from a particle in 
circular orbits around a rotating black hole [9- 10]. 

Strictly speaking, the technique of constructing solutions as series of special functions 
is not new, but extremely old. The similar expansion has been considered in the case 
of Heun's equation[ll-12]. It may not be difficult to see the correspondense between 
the Teukolsky equations for the Kerr black hole and the Heun's equation. That is, the 
Teukolsky equation with two regular and one irregular singularities will be given as a 
confluent case of the Heun's equation which has four definite singularities as we can see 
from Eq.(2.2) in Ref.6 and page 27 of Ref.l2. 

It is known that the separability holds true for the Kerr-de Sitter geometries and even 
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for the Kerr-Newman-de Sitter geometries except for electormagnetic and gravitational 
fields. Therefore, it may be meaningful to consider whether the Teukolsky equations for 
these geometries can be transformed into the Heun's equations, although these Teukolsky 
equations have five regular singularities. A usuful suggestion for this problem can be seen 
in the radial part of Teukolsky equation for the de Sitter geometries [13]. We showed that 
although this equation has four regular singularities, the one of the singularities can be 
eliminated by a suitable change of variables and by the redefinition of the radial function, 
so that the solutions of the equation are given by the hypergeometric function. From this 
experience, we expect that the similar non-trivial transformations may exist even for the 
Kerr(-Newman)-de Sitter geometries. 

Our main aim of this short letter is to show that the Teukolsky equations for Kerr- 
Newman-de Sitter geometries can be transformed into the Heun's equations. Here, we 
stress again that our discussions for the Kerr-de Sitter geometries appliy to massless parti- 
cles, but those for the Kerr-Newman-de Sitter geometries do apply to massless particles ex- 
cept for electromagnetic and gravitational fields. Then, we can use very old technique[12] 
for the analysis of the solutions. The previous solutions for the Kerr geometries [6- 7] and 
de Sitter geometries[13] can be obtained by confluent limits of the Heun's functions. 

In section 2, we provide a short review of the Teukolsky equations for the Kerr- 
Newman-de Sitter geometires. In section 3, we show that both angular and radial equa- 
tions can be transformed into the Heun's equations. In section 4, we analize the solutions 
of the angular and radial equations by using the results of the Heun's equation and show 
that the previously obtained solutions for the Kerr geometries and de Sitter geometries 
can be obtained by the confluent hmits of the soluitions for Kerr (-Newman) -de Sitter 
geometries. The last section is devoted to some discussions and further possible apphca- 
tions. 

2 Teukolsky equation for the Kerr(-Newman)-de Sit- 
ter geometry 

We consider the Teukolsky equations for the Kerr-Newman-de Sitter geometries. Since the 
electromagnetic field and the gravitational field couple through electric charge of black 
hole in this geometry and the equation is not separable so far, we deal with massless 
fields with any spins except for these fields. In the limit of null electric charge, Q ^ 0, 
the Kerr-de Sitter geometries are realized where the Teukolsky equation is valid for all 
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massless fields. In the Boyer-Lindquist coordinates the Kerr-Newman-de Sitter metric 
has the form, 



2 , de^\ Agsin^^ 



+7T-^V^(^^ - «sin2^d<^)2, (2.1) 



where 



Ar = (r^ + a^) (1 - - 2Mr + Q'^ = --^(r - r+)(r - r_)(r - r')(r - r'_), 

Ag = l + Q;cos^, Q; = — 

p = r + ia cos ^, — pp* . 

(2.2) 

Here A is the cosmological constant, M is the mass of the black hole, aM its angular 
momentum and Q its charge. The electromagnetic field due to the charge of the black 
hole is given by 

A^dx^ = . A dt - asin^ dd^). (2.3) 

(1 + ayp'^ 

In particular we adopt the following vectors as the null tetrad, 

^((l + a)(r2 + a2),-A„0,a(l + a)), (2.4) 



V sm6' / 



Assuming that the time and azimuthal dependence of the fields has the form e 
the tetrad components of the derivative and the electromagnetic field are 

v2p V2p* (2.5) 

A/ 2p2' 



where 



A^ A^ ' 



n 

r 

'~'n 



Or + ^ — 7 — V n- 



A, 



A, 



(2.6) 



^ (l + a)H dei^^esim 
Oe H h n- 



A, 



'Ag sin6' 



'l + a)H SefVA^sini 

■ + n- 



A. 



^Ae 806* 



with = cufr^ + a^) — am and if = —acu sin ^ H — : 



m 



sin 6' 



Using the Newman-Penrose formahsm it is known that perturbation equations in the 
Kerr-de Sitter geometry are separable for massless spin 0, |, 1, | and 2 fields. Similarly 
in the Kerr-Newman-de Sitter space those for spin 0, | fields are also separable. The 
separated equations for fields with spin s and charge e are given by 



I AeC\_J AeCs - 2(1 + a){2s - l)aujcos9 



-2a{s - l)(2s - 1) cos^ ^ + A Ss{9) = 0, (2.7) 



ArViVl + 2(1 + a)(2s - l)iujr - ^(s - l)(2s - 1) 

a 

-2{l+a)eQKr + iseQrdrAr + e^Q^r^ ^. ^ ' 
H 2iseQ — A 



Rs{r) = 0. (2. 



3 Transformation of Teukolsky equation to Heun's 
equation 

In this section we show that the Teukolsky equations ( p.TD and (|2.8| ) can be transformed 
to the Heun's equations by factoring out a single regular singularity. 

3.1 Angular Teukolsky equation 

From eq. (|2.7| ), the angular Teukolsky equation become 

I -^(1 + ax^){l - x^)4- + A - s(l - a) + Il±^^2 _ ^^^2 
\ dx dx ^ a 



+ 



1 + a 



1 + ax^ 

(1 + afm'^ 
(l + ax2)(l-x2) 



2s{am — (1 + a)0^ ~ 



:i + «) 



a 



[1 + a)(s^ + 2smx) 
1 — x"^ 



-^2 _2m(l + a)^ + s^ 



>S{x) = 0, 



(3.1) 



where x = cos 6 and ^ = au. This equation has five regular singularities at ±1, and 
00. We also note that the angular equation has no dependence on M and Q. By choosing 
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the variable z 



then eq.(p.l|) takes the following form, 



1 - ^ 



+ 



1 1 1 

+ T + 



z z — 1 z — z^ z — z„ 



x+1 



d 

dz 



f m — sVl f m + s 

T2 



V 2 



+ 



z^ \ 2 
4a 



1 



2 /, . 2 



^1 + ^/am + is' 



s / 1 — u/o' 
+ — 



2ims^/a 



1 + ^^/a 

A - s(l - a) - 2Qi + 2(1 + a)(m + s)^ 



+ 



4a 



2 V^'^(l-^v^)' 



/l+i-\/a\^ 2%ms^f(x 



A - s(l - a) - 2a + 2(1 + a)(m - s)^ 



(1 - z7a)2 



z- 1 



— m 



2ioL^f(x 2^^fa{\ — ol) 



a 



(1 + a) 



(1 + a) 



ms- 



1 + a 



+ 



- s(l - a) + 2) 



{2i^/a.m + (a — l)s)C 



1 + a 

^here z. = -lill^ and 



(l + a)2 
■^i\/a 1 



2; — l + a^; — 
i(l + 



>^(z) = 0, 



-2 - -2o 



(3.2) 



It turns out that the singularity at 



Ay/^ ^ 2v^ 

z = Zoo which corresponds to a; = 00 can be factored out by the transformation, 



S{z) = z^-{z - l)^-{z - z,r-{z - z^)f{z) 



(3.3) 



Tfi — s\ TTl ~\- s\ % I \ ~\~ Ot \ 

where Ai = — - — , A2 = — ;:; — and A3 = ±- ( — — y/am — is . Now f{z) sat- 



isfies the equation 



d^ 

dz^ 



+ 



2^1 + 1 , 2^2 + 1 ^ 2A3 + I 



z-1 



z - z.. 



"2 V 



d ^ p^p^z + u 



dz z{z — l){z — Zg) 



f{z) = 0, (3.4) 



where 



p± = Ai + A2 + A3±Al + l, 



(3.5) 



—I 



u = —^< A - s(l - a) - 2i^ + 2(1 + a)(m + s)i - (1 + i^f{2AiA2 + A + A2) 
4va \ 



4iv^(2AiA3 + Ai + A-i 



m 



4a + (1 + iy/af 



+ —(1 — iyfaf + 2ims\fa{\ + i\fa) 



(3.6) 



Equation (3.4) is called the Heun's equation which has four regular singularities. The f{z) 
is determined by requiring non-singular behaviors at 2; = and 1. We can take either one 
of signs of A3 to find the solution S{z). 



3.2 Radial Teukolsky equation 



From eq.(2^), the radial Teukolsky equation is explicitly written by 



2a 1 

+ 4is(l + apr - —{s + l)(2s + 1^ + 2s(l - a) - 2iseQ - A \R = 0. (3.7) 



This equation has five regular singularities at r±, r'^ and 00 which are assigned such that 
r-t M ± -\/M2 — — Q"^ = r\ and ^^e limit a (A — > 0). By using 

the new variable 



eq (p.7|) becomes an equation which has regular singularities at 0, 1, 2;^, ^oo and cx), 









— r_ \ 




— r_ j 







r_i — r 



r+ — r_ 

Again we can factor out the singularity a.t z = Zoo by the transformation as 
R{z) = z^^{z - lf\z - Zrf^z - z^f^^^g{z), 



(3. 



where 



Si 



5, 



— S ± 2 



2(l + a)a2 (^cj(r2 + 



a2) _ am - ^ 



a(r!|_ — r_)(rL — r_)(r_(_ — r_) 
2(1 + «)a2 (^cj(r^ + 



2\ eyr4 

a — am 



a(r',_ — r_|_)(r'_ — r_,_)(r_ — r_,_) 
2(1 + a)c? [uj{r\^ a') - am - 



IS 



— IS 



a{r^ — r'^){r'_ — r^)(r+ — r+) 



2S 



Then, g{z) satisfies the Heun's equation as 



dz"" 



+ 



25i + s + 1 ^ 2^2 + s + 1 



2^3 + 5 + 1 



z-l 



d 

Z — Zr i dz 

aj^a^z + V 

z{z — 1){Z — Zr) 



g{z) = 0. 



(3.9) 



6 



where 



-s ±i 



2(1 + a)a'^ ( u;{r'J + a^) - am 



l+a 



Q;(r+ — r^)(r_ — r'_){r'_^ — r'_) 
2a*{l + af (r+ - r^)^(r+ - r'_Y{r_ - r'_)(r^ - r'_) 



IS 



r+ — r_ 



I — a;^ri(r+r_ — 2r+r^ + r_r^) + 2auj{auj — m)r-{r-i^r'j^ — r^_) 



'{acu — m)^(2r_ — r+ — r'^) 



eQ 



+ j-j- — [ a;r_ (r+r_ + r_ — 3r+r^ + r_r^) — a(aa; — m)(r+r_ — 3r_ + r+r^^_ + r_r^) 



+ 



eQ 
l + a 



r_ — r 



+ 



2isa2(l + a) 



uj(r_r' + a^) — am - 



a 



(r+ — r_)(r^ — r'_)(r_ — r^_) 



+ (s + l)(2s + l) 



2r 



/ 2 



(r+ — r_ ) (r+ — r'_ ) 

2Bi{ZrB2 + B3) - (S + 1) [(1 + Zr)Bi + ZrBi + S3] 

^ [-A - 2iseQ + 2s(l - a)] . 



a(r+ — r_)(r+ — r'_) 
Here "D is the discriminant of = 0, 



(3.10) 



V 



(r+ — r_)^(r+ — ri)^(r+ — r'_)^(r_ — r^)^(r_ — r^_)^(rl — r'_y 



16a 
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{1-af - {1 - a){a'' + Q^) 



+- 



a 



-27M^ + 36(1 - a)M\a' + Q'') - 8(1 - afia"" + Q'^f 



2\3 



(a^ + Q') 



The sign ambiguity in B2 or i?3 arc related to the boundary condition at the horizon 
or at the de Sitter horizon, respectively. We can take either one of signs of B^. 



4 Solution of Teukolsky equation 

In the former section, we have shown that both the angular and radial Teukolsky equations 
can be transformed into Heun's equations by factoring out the singularity at Zoq. Thus, 
both equations can be solved in a similar way by taking into account the boundary 
conditions. First we briefly explain the formal prescription solving the Heun's equation by 
a series expansion of the hypergeometric functions and then consider each case including 
the boundary conditions. 
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4.1 Solution of Heun's equation 



dz z[z — l)[z — an) I 



The standard form of the Heun's equation with regular singularities at 0, 1, an and oo is 
given by [12] 

7 6 eld aBz — q 

- + r + 

z z — 1 z — an 

where parameters should satisfy the condition, = j+S+e and we assume {anl > 1- 

Solutions which are analytic in some domain including two singularities are called Heun's 
functions. We particularly explain a formal construction of a Heun's function relative to 
the points 0, 1 by means of a series of hypergeometric functions. 

We first construct a solution valid in the neighborhood of 0. The solution is analytic 
in the interior of the ellipse with foci at 0, 1 if the condition for augmented convergence 
presented bellow is satisfied. To the end we write yu{z) in the following form, 

+ 00 

VA^) = E cX+n(^), (4.2) 

n=— oo 

u^{z) = F{-u,iy + uj;r,z) , (4.3) 

where uj = j + S — l = a + P — e. uj is usually used for the definition of Heun's equation 
and should not be confused with the frequency which we use in this paper frequently. 
By using the following recurrence relations[12], 

,, (z/ + cj)(z/ + 7) ,, 2z/(z/ + cc;) +7(cj - 1) ,, 
zuAz) — —- -u^+Az) + zuAz) 

"(" + ^-^1 (4.4) 



{2v + uj){2u + UJ - I) 

d uiv -\- uj){y -\- ^) viu -\- uj){'-y — 5) 

z{z — l)—Uy{z) = —7- -vii^+ii-s) + 7- -^7- -^uJz) 

^ ' dz ^' (2i/ + a;)(2i/ + a;+ 1) (2i/ + u; + l)(2i/ + u; - 1) ^' 

Uu-i{z), (4.5) 



{2v + uj){2v + UJ - I) 

it is shown that yi,{z) becomes a solution if coefficients c„'s satisfy the following three 
term recurrence relations. 



where 



«+i + + 7:<_, = 0, (4.6) 



^ _ {v + n+l){v + n + u-a + l){v + n + u- P + l){v + n + 5) 

~ (2z/ + 2n + c; + 2)(2i/ + 2n + a; + l) ' ^ ' ' 

P'^ = 7 --aH(iy + n)(iy + n + uj)-q, (4.8) 

^" {2u + 2n + u;+l){2u + 2n + u-l) ^ ) \ ) 

^ _ {v + n + a- l)(z/ + n + /3 - l)(t/ + n + 7 - l)(t/ + n + a; - 1) 

~ (2i/ + 2n + a;-2)(2i/ + 2n + a;-l) ' ^ ' ' 



and 

= e(z/ + n)(z/ + n + cj)(7-(5) 

+ [ (z/ + n){u + n + uj) + ap][2{i^ + n){u + n + cu) + -f{uj - 1) ]. (4.10) 

By defining the continued fractions 



RM = -j^, LM = 3;^, (4.11) 

tliey satisfy 



-'n—l '-'n+1 



Now we can evaluate the coefficients by using either series Rn{i^) or L„(i/) with ap- 
propriate initial data. The convergence of the series requires a following transcendental 
equation, 

i?„(z/)L„_i(z/) = 1. (4.13) 

This equation determines either the eigenvalue A for the angular equation, where the 
characteristic exponent u (the renormalized angular momentum) can be taken zero, or u 
for the radial cases, respectively. 

It will be worthwhile to note that if z/ is a solution of eq.( [4.13|) , then — i/ — a; is also 



the solution. This can be shown as foUowings: First, we notice the algebraic relations, 
a-j^-'^ = ryu^ pZn~'^ = Pn and jZn~'^ = aj;. From these relations, we find that cZn~'^'s 
satisfy the same recurrence relations as these for cj^'s. If we choose Cg = Cq "^j we have 
cj; = cZT'^ so that 

Rn{-iy - uj)Ln-i{-y - uj) = i?„(z/)L„_i(z/) = 1. (4.14) 
From eq.(4.3), we see that y-y-u}{z) agrees with yu{z) by choosing Cq = Cq"""^. 

4.2 Solution of the angular equation 

By comparing Heun's equation (4.1) with eq.(3.4), parameters to define Heun's equation, 
a, P, 7, 6, e, u, q and are given by 

a = P+ , /3 = p~ , 7 = 2Ai + 1 , 

5 = 2^2 + 1 , e = 2A3 + 1 , uj = 2{Ai + A2) + 1 , (4-15) 
q = -u , an = Zs , 
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We require that f{z) is a regular function at z = 0, 1, which is satisfied by taking u = 0. 
Then the solution is given in the form of the series of the Jacobi polynomials extending 
from n = to oo, similarly to the Kerr geometry case, 



n=0 

Ur.{z) = F{-n, n + u;r,z) = (_)» r(2n + ^)n! ^(,_,,,_,) _ 

1 (n + 7) 



(4.16) 
(4.17) 



The coefficients determined by solving the equation for Ln{i^) in eq. (|4.6| ) with 

L_i(z/) = 0, where gives the overall normalization. The eigenvalue A can be obtained 
by transcendental equation ([4.13|) with the condition _RJ^ = L^^ = 0. We give A up to 
0(^^) and the coefficient of each power of ^ up to O(a^) as 



A 



/(/ + 1) 



+a 



-1(1 + 1) 



s + 2m^ + 



uJTT) 



l^H{l) + (/ + ifHil + 1) 



+ <^ -2m 1 + 



-2ma 



1(1 + 1] 

o2 



1 + 



/(/ + 1) \ 

+ i^H{l + l)-H{l) + a 



(i-m + i), 



H{1) 



1 + 



-lH\l) + (l + l)H\l + l)- 
H{l + l) 



1(1 + 2) 

H(l + 1) - H{1) + 2((/ + ifHil + 1) - t^H 
H{l)H{l + l) 



Hil + 1) 



+ Qm^s'^(^ 
+ 4m^s^( 



1(1 + 1] 
H{1) 



1(1 + 1)^(1 + 2) (l-l)P(l + l) 
H(l + 1) H(l) 



where 



/2(/ + l)2(/ + 2)2 (/- l)2/2(/ + l)2 

^^^y 2(P-m')(P-sr 



(4.18) 



(4.19) 



(2/-l)P(2/ + l) ' 

and I is the angular momentum which takes an integer or half integer number satisfying 
/ > max(|m|, |s|). From the result it seems that A — s(l — a) is an even function of s. 
Indeed this holds true because eq.( [4.13|) which determines A contains an, Pn and 7„ in 
the forms of a„_i7„ and + 4^(A — s(l — a)), which are invariant under s —s, 
respectively. 

In the rest of this section, we consider two limits, the Kerr-Newman limit and the 
Reissner-Nordstrom-de Sitter limit. It may be worthwhile to note that the angular Teukol- 
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sky equation (3.1) does not contain mass and charge parameters of the black hole so that 
the above limits are only meaningful ones. 
The Kerr-Newman limit: A = (3/a^)a 
In the limit, parameters behave as 



p+ = Ai + A2±s + l + 0{^), p_ 
i 



e + 0(l), A, 



=e + o(i), 



u 



i 



a' ' ' 1\foL 

A - 2A1A2 -Ai-A2 + 2{m + sT (2Ai + 1)) ^ s 

i 



(4.20) 



+ 0(1), 



+ 0(1), 



then parameters in the recurrence relation (|4.7|) — (f4.9|) become by omitting the superscript 



u as 



ar 



i {n + l){n + Ai + A2T s + l){n + 2A2 + 1) 
■7« 2(2n + 2Ai + 2A2 + 3)in + Ai + A2 + I) ' 



(4.21) 



+ 



Jn 



\ ^ 2{n + Ai + A2 + l){n + Ai + A2) 
n{n + 2Ai + 2A2 + 1) 



(4.22) 



A - 2A1A2 - Ai- A2 + 2{m + ST (2Ai + 1))^ 



— s 



In ■ 

where 



i (n + + ^2 ± s){n + 2Ai)(n + 2Ai + 2A2) 
^v^^ 2(2n + 2v4i + 2A2-l)(n + Ai + v42) 



(4.23) 



in = n{n + 2Ai + 2A2 + l){Ai- A2) 

+ (Ai + A2 ± s + 1)[ n(n + 2Ai + 2A2 + 1) + (2Ai + l){Ai + A2) 

Thus, the recurrence relations for obtained from those for c„'s in eq. 



by 



dropping the infinite factor —=. This recurrence relation agrees with those in Ref.4 for 

'a 



the Kerr geometries. In this limit, the Heun's equation ( |3.4| ) reduces to a confluent Heun's 
equation which has two regular singularities at z = 0, 1 and an irregular singularity at 00. 
The angular wave function S{z) is given by using f{z) as 



S{z) 



4\/a> 



(4.24) 



The Reissner- Nordstrom- de Sitter limit: a —* 



11 



In this limit, we have to take a — simultaneously, in order to keep K = "ia/o? fixed. 
Then, parameters become the following forms; 

u = ^^^|\[\-2A,A,-A,-A,-'^^^~s'\J , (4.25) 

'''~4^VA' ''°°~2^VA' 
Then, and 7^ are 0(1), while [3^ are 0{l/a) —>■ 00, 

i [3 



0" 

4a V A 



w? — 



n{n + 2Ai + 2^2 + 1) - A + 2A1A2 + Ai + A2 + ^ + s 



(4.26) 



In order that coefficients a„ have definite limits, must be for a specific value of n. 
The eigenvalue A is determined by this condition as 

A = (/-s + l)(/ + s), (4.27) 

where I = n + Ai + A2 > max(|m|, \s\). Therefore the angular wave function has the 
following form, 

xF(-/ + A^ + A2,l + Ai + A2 + 1; 2Ai + l]z) . (4.29) 



In this limit, the Heun's equation ( p.4| ) reduces to a hypergeometric equation. These 
results coincide with that of [14]. 

In the former subsection, we comment that if z/ is a solution of eq.(4.13), then —v — uj 
is a solution. For the angular solution, we chose z/ = and we know that a;=integer so 
that this transformation of the exponent does not lead any new solution. 

4.3 Solution of the radial equation 



By comparing eq. ( |4.1| ) with eq. (|3.9|) , we find parameters in Heun's equation, a, /3, 7, 5, e, q 
and an is given by 

a = (T+, P = (T_, 

7 = 25i + s + l, 6 = 2B2 + s + l, e = 2^3 + 5 + 1, (4-30) 
q = -V, an = Zr. 
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The solution convergent in the eUipse with foci at ^ = 0, 1 which correspond to r = r_, r_| 
respectively is given by 



Uu{z) 



+00 



F{-u, V + 2(5i + 52 + s) + 1; 25i + s + 1; z). 



(4.31) 
(4.32) 



where to' = 7 + i5 — 1 = 2{B\ + ,82 + s) + 1- The expansion coefficients fej^'s are determined 
by the recurrence relation (|4.6| ) in which c„'s are replaced with foj^'s. The radial solution 
is expressed by a series where n runs from —00 to +00. This is because the eigenvalue A 
is fixed from the angular solution and thus the characteristic exponent (the renormalized 
momentum) v is determined by the transcendental equation ( |4.13| ) which guarantees the 
convergence of series. Remind that if gy{z) is a solution then g^v-^{z) is a solution too. 
the Kerr-Newman limit: A = (3/a^)a — >■ 
In this limit, r± and are 



r± 



1 + 



2a^{rl-M) 



1 ± — V" H —r> OL 



then parameters are 

5i = ^[-sTi(-e + f + is)]+0(v^) = fi° + 0(v^), 

B2 = i[-s±2(e + f-2s)]+0(v^) = 5° + 0(v^), 

„ iauj 1 , , 

= T-^ + -[-s±«(-e-«s)]+0(v^), 

C7+ = 5° + 5° + s + l±i[-e-is] +0(a), 

o = ^!^ + 50 + 50 + , + i + o(v^), 



+ 0{a), 



(4.33) 
(4.34) 



(4.35) 



2^/a{r^ — r 
-iek 



0- {-A + 2s- 2iseQ + 25°5° + (s + 1) (b^ + 5° 



+0(v^) 
a 



s ± (25° + s + 1)] + 2e2 - t^k + ege(K - 3) - eg - ^ (e^ + - 2isi) | 
+ 0(1), Zoo = ^^-^ oT + 0(l)' (4-36) 



2ya(r0 - rO ) 



^J^(^r% — r 



where e = 2Mu - eQ = e - eQ, e = 2Mu;, k = Jl - <? = """t? and f 



e - q 



The parameter u in the above expressions is used for the frequency. The singularity at 
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CO becomes irregular because of oo. Thus the radial equation (|3.9| ) becomes to a 

confluent Heun's equation as the angular one does. In the Kerr limit (A ^ oo and Q = 0), 



coefficients of the recurrence relation 
a 



coincide with those in Ref.6, if we divide them 



the de Sitter limit: M, Q and a — > 

We rearrange the locations of the singularities so that z = 0,1, z'^, and oo correspond 
to r = r^, r_, r_|_ and r'_ respectively; 

(r_ — r'_) (r — r',) 



(r_ — r+) (r — r'_) ' 



(4.37) 



thus zL 



fr_ — r' ) fr I — r'. 



The effect of this rearrangement is only to replace r_ , r_|_ , r'_^_ 



and r'_ with r', , r_, r+ and r'_, respectively, in the radial equation. In this limit we find 



r± — >• ±ia. 



(4.38) 



and 

B2 



-uj + is 



- [-S ±{m + s)] , Bs 



- [-S ± (-m + s) 



-s ±i \ ^' ~~ "^-^ 



-A + 2s - (s + l)(2s + 1) - 2Bi{B2 + ^3) - (s + l)(25i + B2 + B3) + is 
1 z = —1 



(4.39) 



-Co", 



The resulting radial equation is the following with the three singularities 

r d-" 



dz^ 



+ 



2Ei + g + 1 ^ 2(^1 + B2 + S + I) 



d a^a_ + f 
dz z{z — 1)^ 



^(z) = 0. (4.40) 



z z — 1 

By using A = (/ - s + 1)(/ + s) ( ^1271) and setting ^(2) = (2 - l)"'^'"^'"'"^+^'+5)^j^g(^) 
this equation reduces to a hypergeometric equation thus we obtain 



R{z) 



5'dcsl2;j 



[z + l)2^+iF(/ + 1 - s, / + 1 + z J-cu; 1 - s + ^y-cu; z).(4.41) 

—UJ because of a 



This radial function coincides with that in Ref.l3 except for uo 
difference between our choice of the null tetrad and theirs. 
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5 Conclusions and discussions 



In this note, we have shown that equations for perturbations of the Kerr-de Sitter black 
hole for massless particles with any spins are reduced to the Heun's equation and thus we 
can obtain solutions in the form of series of hypergeometric functions. This method has 
been extended to the cases of Kerr-Newman-de Sitter geometries except for photons and 
gravitons for which the technique to separate angular and radial parts are not succeeded 
so far. We have explicitly constructed solutions for angular and radial equations and 
examined some limiting cases. We have shown that these limits reproduce the former 
known results. 

The radial solutions we obtained in the text are valid inside a ellipse with foci at 
z = 0, 1 and are not valid at de Sitter horizon. In order to apply for some physical 
situations, we have to construct solutions which are valid in the entire regions of r and 
satisfy some specific boundary condition. For this, we should obtain solutions valid around 
the de Sitter horizon and match solutions with different convergence regions in the region 
where both solutions are convergent. This procedure was made for the Kerr black hole 
case in Ref.6,7. Also, it may be interesting to see the Teukolsky-Starobinsky identities 
analytically. We will answer these questions in the future publication. 
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